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The joint ES A/NASA LISA mission consists in three spacecraft on heliocentric orbits, flying in a 
triangular formation of 5 Mkm each side, linked by infrared optical beams. The aim of the mission 
is to detect gravitational waves in a low frequency band. For properly processing the scientific 
data, the propagation delays between spacecraft must be accurately known. We thus analyse the 
propagation of light between spacecraft in order to systematically derive the relativistic effects due 
to the static curvature of the Schwarzschild space-time in which the spacecraft are orbiting with 
time-varying light-distances. In particular, our analysis allows to evaluate rigorously the Sagnac 
effect, and the gravitational (Einstein) redshift. 
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I. INTRODUCTION 

The three spacecraft forming the LISA constellation 
have a complex motion (see for instance) allowing the 
triangular constellation to remain approximately rigid 
during the annual revolution. There is simultaneously 
a rotation of the triangle around its center, and the or- 
bital motion of its center. A number of papers have con- 
sidered the effect of the rotation on light propagation. 
Accurate knowledge of light propagation delays are very 
important for implementing the TDI (Time-Delay Inter- 
ferometry) technique (see 0) which is mandatory for 
eliminating laser phase and similar noises. In particular, 
in a simulation code, it is necessary to generate time de- 
lays as realistic as possible. This is why we have carried 
out a full general relativistic treatment of light propaga- 
tion in the gravitational potential of the Sun. The back- 
ground metric created by the Sun is responsible for both 
the orbits, and the fine structure of time delays. Even 
if the orbit of any spacecraft is taken as purely Keple- 
rian, the time interval between the emission of a photon 
off spacecraft A and its detection in spacecraft B must 
be thoroughly investigated. In particular, the motion of 
the target B during the interception must be evaluated 
with an accuracy level consistent with that of the null 
geodesic followed by the photon. We therefore propose 
an expansion of both the null geodesic equation from A 
to B and the motion of the target B, in powers of the 
small dimensionless parameter 

e = rsl2r = GM/{rc^) (1) 

where rg = 2GM / (? is the Schwarzschild radius of the 
Sun (of mass M) and r the coordinate distance of the 
spacecraft from the Sun. We intend to derive the time 
delay between spacecraft and the global frequency shift 



(including classical Doppler due to relative motion or or- 
bital shift plus Einstein shift). The first quantity is es- 
sential for TDI, the second seems important for technol- 
ogy, because large frequency shifts, variable in time over 
a year induce requirements on the ultrastable oscillators 
used for compensation of orbital frequency shifts. We 
propose formulas for computing the time delays neces- 
sary for a realistic numerical model and we show that in 
the global frequency shift, the overall general relativistic 
contribution (first order terms) is much smaller than a 
naive estimate would suggest. 

The method we present is based on a matching of 
geodesic equations of spacecrafts A and B, with the (null) 
geodesic equation of the photon emitted by A and re- 
ceived by B. 

II. COMPUTATION OF (COORDINATE-)TIME 
TRANSFER IN THE LITERATURE 

Different methods have been considered for computing 
the flight time of photons and the gravitational frequency 
shift between two space-time points x\ = {tA,^A{tA)) 
(emitter) and x% = {tB,^B{tB)) (receiver). An alterna- 
tive to integrating the geodesic equations of motion is to 
use the world function Q (x^, x^) developed by Ruse 
and Synge Q, defined as the squared geodesic distance 
between two space-time points with r2(a;^,a;") = be- 
ing the equation of the light cone at x^. The time delay 
corresponds to the future cone equation, and is given 

(PAT) 

by T {tA,XA{tA),XB{tB)) = ts ~ tA ^ Rab/c- Vl 
(tyi, x^(i:^), XB(i_B)) /(c Rab), the subscript PN mean- 
ing Post-Newtonian contributions, Rab = ^B{tB) — 
:x.AitA) and Rab = \'Rab\- 

Linet et al (2002) provide the world function. 
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V, {tA,XA{tA),^B{tB)), up to order 1/c^ in the full 
Nordtvedt-Will Parametrized PN (PPN) formalism (10 
parameters characterizing alternative theories of grav- 
itation, including the usual PN parameters (3 and 7). 
Their results yield an expression for the coordinate time 
transfer, t {tA,^A{tA),^B{t-B)), at order 1/c^. The au- 
thors also provide r (iyi, Xyi(tyi), XB(iA)) which contains 
explicitly the motion of the receiver and the correspond- 
ing Sagnac effects of orders l/c^, l/C*. They ap- 
plied the earlier expression to an isolated, axisymmetric 
rotating body, assuming a stationary gravitational field 
and a constant velocity of the central gravitational body 
with respect to the universe rest frame. Their systematic 
procedure to compute multipole moment contributions 
in (xA(iA),XB(tB)) and r (xA(iA), XB(tB)) was then 
particularized to obtain the explicit contributions of the 
mass monopole, mass quadrupole moment and of the in- 
trinsic angular momentum of the rotating gravitational 
body. Restricting to fully conservative metric theories 
of gravitation, without preferred location effects (that is 
all PPN parameters vanish except f3 and 7), they de- 
termined the fractional frequency shift up to order 
The authors then performed numerical estimates of the 
frequency shifts in the gravitational field of the Earth. 
They assumed A on board the International Space Sta- 
tion (ISS), orbiting at an altitude of 400km, and B in 
a terrestrial station; as this is the case of ESA's Atomic 
Clock Ensemble in Space (ACES) mission, planned for 
2009-2010, aiming at an accuracy of 10^^^ in fractional 
frequency. The formula the authors give yields all the 
gravitational corrections to frequency shifts up to 10~^^ 
in the vicinity of the Earth. 

The above cited paper in fact generalizes the work of 
Blanchet et al (2001) [3, based on the geodesic approach, 
who provided the time transfer and frequency shifts up 
to order l/c"^ in the setting of General Relativity (where 
/3 = 7 = 1 are the only nonvanishing PPN parameters), 
for a monopole mass, without intrinsic spin. Paper [3 
considered the two-way laser link in addition to the one- 
way laser link presented in paper However, for the 
LISA mission, only single links arc pertinent since laser 
light is not reflected from one station to another, LISA 
consisting in three pairs of laser links. 

Le Poncin-Lafitte et al (2004) Q present a generic 
procedure based on the world function to provide, in 
an iterative way, the propagation time of a photon be- 
tween space-time events Xa and x% at the nth post- 
Minkowskian order (developpement in powers of G) from 
the contributions of orders p = n — 1, n — 2,..., 1, 0, 
(p) 

i.e. of il (i^, x^(tyi), XB(iB)). The authors then apply 
their formalism to a static spherically symmetric space- 
time within the second post-Minkowskian approxima- 
tion, with the usual post-Newtonian parameters /3, 7 plus 
parameter S for the G^-term in gij (normalized so that 
in General Relativity 6=1). Within those assumptions, 
they obtain the expressions of the world function and 
time transfer in the case of a simple monopolar and non- 
rotating body, up to order (containing order 1/c^). 



The motion of the stations (and corresponding Sagnac 
effects) is not explicitly considered there. 

We are interested in the time variations of the inter- 
spacecraft propagation delays over a year for spacecraft 
having very special orbits, this is the reason why we de- 
velop a special approach, on principles analogous to the 
preceding works, but in a form more suitable to our pur- 
pose. 

III. POST-NEWTONIAN THEORY 
A. Notations 

Coordinates {ct, x, y, z) are denoted by x". We use the 
notation r = \/~xFx} ^ where latin indices take only the 
spatial values 1, 2 and 3. Then, r cannot be understood 
as a physical distance in all the steps of the developments. 
It should be referred to as a coordinate distance. 

Four-velocities are written u" = da;"/dr, where r is the 
proper time, given by c^rfr^ = —ds^, since the signature 
chosen for the metric is (— , +, +). We frequently use 
the three- velocity = dx'^/dt. This quantity cannot be 
understood as a physical velocity in all the steps of the 
developments. When a quantity Q is evaluated up to 
order p (integer or half-integer) in e, it will be written 

(^p) (0) (1/2) (p) 

Q =Q + Q +...+ Q 

(0 . 

where Q is the contribution of order I. 

The notation U represents the three numbers , and 
U.V represents the sum WV^ = U^U^ + U^U^ + U^U^ . 
The notation = U.U will also be used, rjap = 
(—1, -1-1, -1-1, -f 1) is the Minkowski metric. 

In this paper, we derive the (relative) frequency shift 
of a photon, linking an emitter spacecraft A to a receiver 
spacecraft B, up to order 3/2 in e. 

B. Generalities on the metric 

We consider a metric of the generic static form 

ds^ = -I{e)c^dt^ + J{e)6ijdx'dx^ (2) 

where / and J are two functions depending on the theory 
(e.g. General Relativity, Scalar- Tensor theory,...). Both 
are assumed expansions in integer powers of e. The lead- 
ing order in f'/c, is 1/2 (Kepler). A photon worldline, 
described by the wave vector A:" = dx'^/dX, where A is 
an (arbitrary) affine parameter, obeys the isotropic con- 
dition 

kak" = (3) 
and the null geodesic equation 
dk°' 

— -i-r^,fc^r = (4) 
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or, as well 



dk 1 



Since the considered metric is stationary, fcp = ~I{f)k^ 
is a constant, according to ©. 

A free falling observer of coordinates a;" and four- 
velocity u" can be described by the geodesic equation 



D. Light propagation in curved space-time 

(5) 

The metric tensor is of the form 

gap = riap + 2gQ e S^p (9) 

where go = 1 and qi — "f. We have as well 

5"'' = Vo.p-iqaeS"" (10) 

The motion of the photon is described by the geodesic 
(6) equation 10} or Q and by the null condition (Q. 



C. Energy of a photon 



Order 



Consider a photon of four-wave vector fc" propagating 
in the background metric Q. The general expression of 
its energy measured by an observer of four-velocity 
reads: 



E = -gapk^u^ = -fcow° - J(e)k.u. 



(7) 



An expansion of E up to order 3/2 requires a knowledge 
of vP /c and /c at the same level. On the other hand, 

(0)' 

u'/c being a term of order 1/2 (i.e. u = 0), we need 
A;' and J up to order 1 only, to evaluate E. Since the 
metric components are integer powers of e, the functions 
J(e) and J(e) have to be known up to order 1 only, in all 
the developments we will have to consider. 
Let us compute u*^ — c.dt/dr. We have 

c^dr^ = I{e)c^dt^ - J{e)dx'dx' = (^^ - 

Let us consider metrics of the form 

/(e) = l-2e + 0(e2) 

J(e) = 1 + 276 + 0(e2) 



J- 



c'dt' 



where we recognize the isotropic form of the 
Schwarzschild metric in the case where the Post- 
Newtonian parameter 7 is unity. We thus have at the 
same order 

2e- - 



dT^ 



1 



and consequently 



c 



dt^ 



We notice that proper and coordinate times differ by 
terms in e (Einstein shift) and in w^/c^ (Lorentz time 
dilation). Now we compute 
,,0 



At this order, space-time is flat and the spacecrafts 
are at rest (no velocities, which are 1/2 order terms, and 
no accelerations, which are first order terms), and light 
propagates according to special relativity. The results of 
order are obvious, but since the same notation will be 
used in higher order terms, we present here the detailed 
solution. 

Solving eq.Q), we find, the F being first order terms 



(0)" 

k 



where C" are four integration constants. From the 
isotropic condition (O, 

G/^ /^O/^O r^i r^i 

where Cq, = r^apC^ , so that there is a spatial unit vector 
n* = jC^ along the direction of the (straight) light ray. 
We can parametrize the order photon worldlinc using 
the trivial differential equation 

d X (0) 
~dX~ 

which gives (assuming A = for the emission time) 



(0)" 
X 



and if x\ represents the coordinates of the emitter at 
emission time, we have 

(0)* 

V ^x\ + C'X 

(0)° 

Let us remark that interpreting x as the Minkowski 
time coordinate ct, this could be simply written as 



j dx^ dt dx 

dr dr dt c 
from what we get the energy 
E 



) .vo(.^/^) 



XA + ct.n. 

(0), 



The spatial coordinate distance r (A), which corresponds 
to the "physical instantaneous distance" at this order, is: 



- = -fcn 1 



2c2 



l + (l + 27)e+^ 



k.v 

c 



(0) 



r{XY 



(o)'(o)' 

: X X 



[C°A + n.XA]' + if' 



(11) 



(8) with K"^ — r\ — (ii.xa)^ and r\ = x' 
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2. Order 1 

ko being a constant, we obtain k'^ at order 1 by 

(^1)" (^1) (0) (0) 

'k ^g"' k 0^ r Vo= g™Co - (1 + 2 (12) 

(0) , (0) 

where e = 2rs/ r . Integrating the spatial components 
of the geodesic equation or one finds 



n*-(l + 7) e F*n. x -(7-I) e 



where is defined as 
pi 



(13) 



if2 



and satisfies P.n = 0. Integrating equation 
one finds, in parametric form 



a; (A) = C°A + 2 



(0) (0) 

GM , n. X + r 
— TT- In 



n.XA + TA 



^~t^\x) ^x\ + C"A - (7 + l)^P'tr -rA] 



J, n. X + r 

- 7 - I Wti* In ■ . 

c/ n.XA + 

Ehminating the affine parameter, one obtains the pho- 
ton's trajectory in the exphcit form 

f-*i~i* CM 
x' {t,n^)=x\ + n'ct-{-f + l)^x\t.ri^) (14) 



with 

(0) (0) 

xHt, n^) = P'Vr -ta] + In 

n.XA + TA 

where the dependence in the integrations constants 

is exphcitely written and is given by (|ll|l with C°A 
replaced by ct. 



3. Energy of the photon 
Inserting (|13|l in ©, one finds 



I E _ n.v GM 

CO 7 " " ~ ^ 7^ ^ 2? 



, GM v2 
(2 + ^)(^^ + 2? 



r 



n.v 

c 



+ (7 + 1) 



(0) 

GM P.v n. X 

~^ c (oT' 



(15) 



The first line of the right-hand side (r.h.s.) of (|15|) con- 
tains terms of order (i.e. 1), 1/2 (i.e. — n.v/c) and 1, 
while terms of order 3/2 are in the second and third lines. 

Equation 115|) shows that knowledge of the energy up 
to order 3/2 requires 

- knowledge of n* up to order 1; 

- knowledge of /c up to order 3/2; 

- knowledge of up to order 1/2. 



4-. Evaluating the flight time and the constants 

At t ~ Q, a photon is emitted from x\ by spacecraft A 
with initial position x)^ and velocity v^, and is received 
at ^^(t) by spacecraft B with initial position and 
velocity v^. We want to compute both the flight time t, 
and the constants n* characterizing this photon. 

The equation to solve writes 



x\t,n^)^x^^{t) 



(16) 



where x'^{t,n^) represents the photon's trajectory p4|) 
and x^^{t) the receiver spacecraft's trajectory. Since i = 
1, 2, 3, equation H16(l represents three equations involving 
the four unknown t and n\ The fourth equation is the 
normalization condition 



n^ = n'n' = 1. 
The Taylor development reads 



x],{t) 



tvl 



dt 



(17) 



(18) 



(t=o) 



Since the operator d/dt = v'^di increases by 1/2 the order 
of the quantity on which it is applied, this development 
shows that : 

- knowledge of a:;g(i) up to 1/2 order requires knowl- 
edge of t at order 0; 

- knowledge of v^{t) — dx'^{t)/dt up to order 3/2 re- 
quires knowledge of t at order 1/2. 

Hence, one has to solve equation (|16|l up to order 1/2 for 
t and up to order 1 for n*. 

Let us write ()16|l in the explicit form, using H14|l and 
(EJ), 

x\ + n'ct - (7 + l)—x\t, n')^x\,+ tv], ^ 

C ^ '^B 

(19) 

where the initial acceleration, given by the geodesic equa- 
tion (0 reduces to its newtonian part at order 1. We 
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notice that, due to properties of the metric considered 
(static and diagonal), this is also the expression of the 
initial acceleration up to order 3/2. 

At order 0, equations ((T^ and l(T7|) read 

Xj^+ n c t = 



(if 1 (0)' 

n = n 

2c2 



2 /(O) ^2 

VB-(n .vb)' 



(0) 

GM t 
2r|c 



i (O)'(O) 

x^— n n .xb 



n = 1 



which give the two (obvious) order solutions 



+ (7 + 1) 



GM 

c3 t 



, (O)'(O) 

X- n n .X 



(0) (0)' 

t . n 



(0) 
C t = 



. (1) . ■ . 

The "In" term in t , present in the term, is responsible 
for the well-known Shapiro time-delay effect. 



(0)' , 
n = ± 



(21) 



where the sign corresponds to a photon emitted by 
spacecraft A and received by spacecraft B, while the "— " 
sign corresponds to a photon emitted by spacecraft B and 
received by spacecraft A. This is nothing but the trivial 
euclidean prediction in the case where both emitting and 
receiving spacecrafts are at rest. 

At order 1/2, using the order equations, equations 
(O and ini) read 



(0)'(l/2) (0)(l/2)* (0) vl: 

n t + t n = t — 



(0) (1/2) ^ 

n . n = 0. 



The corresponding terms of order 1/2 



(0) 

(i/2)_(o) n .VB 



(1/2)' 

n 



(22) 



(0)' n .vb 
n 



are the classical motion and aberration corrections. 
The same iterative method provides the next order 

(1)' 

terms. For the energy problem, one only needs n , but 
(1) . . 

the method also gives t , which will be useful in the 
time-delay problem 



(1) 1 

'-2 



(0) 

n .VB 



(0)' 

n 



(1+7)^X ( t,n] 



(0) 

t + 



GMxl (0) 

t 

2r|c 



(23) 



IV. TIME DELAYS 

The time transfer between two satellites A and B is, 
at order 0, 



(-1) (0) 

t = t 



(1/2) 
t 



(1) 
t 



where the contributions t have been calculated in the 
previous section. The order term is the same for both 
links A^B and B^A. The first term to be affected by 
the permutation A ^ B is the 1/2 order term. It is easy 
to see that, for this permutation, one has 



.(0) (0)' (0) (0)\ 

(i,n)-(f,-n) 



leading to 

t (A 



B)^ 



(1/2) ^ 

t (B 



A) 



(0)(0)' 

= t n 



which is responsible for the so-called Sagnac effect. Let 
us notice that in the TDI problem, one generally needs 
the emission time as a function of the reception time (or, 
equivalently, the flight time as a function of the reception 
time), rather than the converse. Accordingly, the solu- 
tion required in simulations generally involves the order 
solution 1201), lETJ with the minus sign. 



V. EINSTEIN FREQUENCY SHIFTS 

Three identical clocks, one aboard each LISA space- 
craft, beat at a common proper frequency. The quantity 
measured on spacecraft B is the (relative) difference be- 
tween the frequency received from spacecraft A, and the 
proper frequency measured on B. For a photon starting 
from A at t = 0, and received at B at time i, this relative 
frequency shift writes 



Z(B^A)(iA = 0) 



EAit = 0) 



1 



(24) 
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since the nominal frequency of the oscillator aboard B is 
identical to the one of the oscillator aboard A, which is 
nothing but the initial energy of the photon emitted from 
A. 

Let us remark that 1 + z as defined in (|24|l corresponds 
to the standard definition of a frequency shift (frequency 
at reception over frequency at emission) in the case where 
the solution (|20I21|I with the + sign is considered. In the 
other case (— sign), 1 + z as defined in H24|l corresponds 
to the inverse ratio (frequency at emission over frequency 
at reception). 



(0)' 

where P = 



n ~ n 



Since l|25|) has no order 



term, one only needs up to order 1 to obtain z up to 
order 3/2. Setting D = E^it = 0)/(cC°) yields 

(0) 

D= 1 



(1/2) n .va 
D = 



(28) 



A. Contribution of the different orders 

The development of the expression for the energy up 
to order 3/2 is given by l|15|) . This expression requires 
the velocity of the spacecraft up to order 3/2 



(3/2)* . ^ dvl. 

V B W="B+i^ + 



tGM-^ - —GM 



"2 di2 



x^Xb.vb 



Substituting the above expression in equation (|15|) . one 
finds, writing /^E = E^iit) - EA{t = 0) 



(1/2) 

AE =AE 



(1) 
AE 



(3/2) 

AE 



with, writing N = AE/{cC° 



(1) 
iV = 



(1/2) 

N 



n — 



(25) 



(26) 



(o)'(o) 1, ^ 

n n .vb - -(wb 



«1) 



GM 



1 1 



GM (0) ^ ^^^^ 



(3/2) _ (0)^ „^V| - VIW\ 

TV = - n — n -2—2 — A A 

c 2c3 



(1) 
D- 



2 \ B 



(O)'(O) 

n n 



GM 



Hence, one finds 



(1/2) 
Z 



(1/2) 

= N 



(1) (1) (l/2)(l/2) 

z^N - N D 



(29) 



(30) 



(3/2) (3/2) 

z = N 



(l)(l/2) 

N D 



(1/2)(1) 

N D 



(1/2) 

N 



'(1/2)' 

D 



(1/2) 



It appears that z is unchanged under the permutation 
A ^ B, as it should be, since the Dopplcr effect only 
depends on relative motion in the framework of special 
relativity. Of course, this is not the case as soon as grav- 
ity is acting, i.e. in first and higher order terms. Indeed, 
one finds 



(1) 

Z (B*- 



-A) (iA = 0)- 



Z (A^B) (.IB = Uj 



2GM 



1 



1 

TA 



GM (0)^(0) 
n t 



The first term originates in the Einstein gravitational 
Doppler effect, while the second is related to the accel- 
erations of spacecrafts in the gravitational field of the 
Sun. 



B. Orders of magnitude for the LISA configuration 



GM (0) vb.xb GM {QfW 



-(2 + 7) 



GM {of 



a;^XB.VB 



_ 3: 



rB TA 



+ (1 + 7) 



GM (o)''(o)^' 
P n 



■-b^b 
tb 



'-A'^A 
''A 



One is tempted to derive an order of magnitude for 
the contributions to the frequency shift. Indeed, one has, 
from 

10-8 

and the difference between the distance of each station 
to the Sun satisfies 

5r L TT 1 
— < - cos - - — 
r r 3 oO 
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where L ^ 5.10^ km is the typical interspacecraft dis- 
tance. An estimation of the terms of orders 1/2, 1 and 
3/2 leads to 



(1/2) 

Z r 



.1/2^^2.10^6 



(3/2) 

z . 



5r 



3/2 ^ 



2.10- 



-10 



2.10 



-14 



It is interesting to point out that, due to the peculiar 
LISA configuration, terms of orders 1/2 and 1 are actu- 
ally considerably smaller than the above estimates. Let 
us evaluate these terms more carefully. 
In the ideal case, where LISA spacecraft form a perfect 
equilateral triangle, velocities of spacecrafts A and B have 

(0) 

the same projection on vector n , such that from H29|l . 



(1/2) 

z ~ Q. 



In fact, the above cancellation is obtained as soon as the 
orbits about the center of mass of the satellite constel- 
lation are coplanar and circular, with a velocity propor- 
tional to the center of mass distance. Since this property 
is not exactly verified by the real configuration, the value 

(1/2) (1/2) 

of z is reduced by a factor ^ L/r, i.e. z ~ 7.10 . 
At this order, the frequency shift is due to the arm length 
variations (flexing). 

Let VAB = VB - VA. From ^ and 

the term of order 1 writes 



(1) 
z - 



(0) 

n .vab 



VAB 



GM ( 1 



^ (31) 



The first term on the rhs is small, ~ 5.10 since it is 

(1/2)^ 

nothing but z 

In the second term, v^g turns out to be three times the 
squared modulus of the relative velocity of a spacecraft 
with respect to LISA's center of mass, i.e. VL/{r^/2>) 
with V the orbital velocity of this center. Hence, 



1/VAB\2 12 

The third term of the rhs of (|31|l is of order 
GM 5r 



2.10 



-10 



while the fourth writes 

(0) 

GM (0) n .XB 

5 C t o 



^^cos0^2.1O- 



10 



where 9 spans between tt/3 and 27r/3. 
The point is that the third and fourth rhs terms of equa- 
tion (|^ . of order {GM / c^){L/r'^) each, nearly cancel. 
Indeed, the euclidean triangular assumption. 
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\ - 2rBL.XB + 
with XA B — f'A.B XA.B and L 



c t n leads to 



1 

''A 



1 

rB 



L.Xb 



-1 



(0) _ 

n .XB 



2rl 



-O 



, (0) _ 
Since ( n .xb 



< 1/4, the sum of the third and fourth 
terms in (|31|l is about 



1 GM 



2 rc2 



6.10^ 



(33) 



Furthermore, it turns out that the residual 
(third-hfourth)-term in of order {GM/c^){L'^/r^), 

nearly cancels with the second term in 1311) . Indeed, at 
first order, each spacecraft obeys Newton's law 

so that their relative motion is described by 



^ = -GM(k-'S^'-^^' 
dt 



with Vab = dL/dt. Since ^(L.vab) = L-c'vab +^^36?^, 
it turns out that the sum of the second, third and fourth 
terms of (|31|l reduces to 



J_d_ 

'2^dt 



L n .VAB 



O 



i3 

^4 



(34) 



The order of magnitude of the 0{L^ /r^) terms neglected 
at the different steps leading to (|34|l is 



GM 



■cos 6* - 2.10 



-13 



But the first term of (|34|l is even smaller, since 



(1/2) 

1 / (0) \ I r Z 



(1) 



3.10" 



(32) Consequently, z appears to be of order 2.10 , i.e. only 

(3/2) 

an order of magnitude larger than z . 

In the next section, we will confirm the partial can- 

(1/2) (1) 

celling of z and z found analytically, thanks to a 
numerical relativistic simulation of LISA's optical links. 
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VI. NUMERICAL RESULTS 

In order to evaluate the orders of magnitude of the 
various terms of the expansions in the time delays and in 
the frequency shifts, we have made a toy model in which 
the spacecraft orbits are classical (Keplerian). At any 
date, the link between spacecraft A and B is assumed 
starting from the Keplerian location of A at current time 
tA and joining B at a place deduced from the keplerian 
position of B at time tA + 1, through a relativistic contin- 
uation. It is known (this has been discussed for instance 
recently in |l|), that the three spacecraft can be kept at 
approximately constant mutual distances, with the con- 
stellation's center of mass on a circular orbit of 1 A.U. 
radius, at the condition that each orbit has a small eccen- 
tricity e, a small orbital inclination ^ on the ecliptic, with 
the semimajor axes of spacecraft 2 and spacecraft 3 orbits 
shifted by 120 and 240 degrees respectively from space- 
craft I's. In fact, if the three spacecraft are required to be 
equivalent (i.e., same orbit eccentricity e and inclination 
angle fj,), at first order in e, the distances between space- 
craft are constant and the angle of the LISA plane with 
respect to the ecliptic, i^, is 60 degrees. However, exact 
numerical investigations show that second order terms in 
e give rise to important variations of the intcrferometric 
arm lenghts ("breathing") and to a large Doppler effect. 
The choice of the orbital parameters (e, ^)influences the 
"breathing" amplitude. A carefuU choice of the orbital 
parameters can reduce the "breathing" amplitude from 
about 120,000 km (for a naive choice of 1^=60 degrees) to 
less than 50,000 km. Such orbital parameters were cho- 
sen in existing simulation codes ^^'^ been 
shown in reference illj that this choice is the optimal 
one. We summarize herafter the results of terms 
of parametrized orbits. The angle v is defined as 

TT 5 

= — I — a 
3 8 

where a = L/2R. L is the nominal interspacecraft dis- 
tance (5 10^ km) and R the radius of the circular LISA 
orbit (1 a.u.), so that a ~ 1/60 is a small parameter. 
This small change of inclination is responsible for a large 
reduction of flexing with respect to the strict 60 degrees 
often quoted in the literature 0|. The inclination /i of 
the orbital planes with respect to the ecliptic is: 



tan /i 



asmz^ 



sin(7r/3) 



a cos ly 



Though the center of mass keeps on a circular orbit, the 
spacecraft stay on slightly elliptical orbits (numbered b = 
1, 2, 3). Their common eccentricity is: 



4a cos i> 

~7f~ 



4q2 
~3~ 



1/2 



In order to parametrize the orbital motions by the time, 
we need the eccentric anomalies ^j,, that are functions of 



the time t, defined via the implicit equations: 

'^b- e sin = ili - Ob 

involving a phase shift 6b = 2(6 — l)7r/3. Owing to the 
smallness of e, the numerical solution is very fast. Then, 
let us set 

Xb = R (cos \I'6 — e) cos /x 
Yb ^ RVT^ sin*b 
Zb = — -R (cos ^Pf, — e) sin/i 

Finally, the positions of the three spacecraft at time t 
are deduced from the preceding ones via rotations in the 
{X, y)-plane: 



Xb — Xb cos 6b — Yb sin 9b 
yb = Xb sin 6b + Yb cos 6b 
Zb = Zb 



(35) 



At this point, we have three Keplerian orbits, with pa- 
rameters such that the distances between test masses are 
constant at flrst order in a, but actually varying by about 
1% (about 50,000 km) peak-to-peak during the year. The 
spacecraft velocities are computed according to: 



Xb = Xb cos 6b — Yb sin 9b 
yb = ^6 sin 6b + Yb cos 6b 
Zb = Zb 



(36) 



with: 




QR 



1 — e cos vj/i. 



— cos ^ sin 'i'b 
\/l — cos ^I'f, 
sin fi sin 



Once the positions and velocities are known, it is easy to 
implement the time delays up to order 1. Fig lll2l3l show 
the contributions of orders 0, 1/2 and 1 respectively. We 
can also evaluate the Sagnac term as the difference be- 
tween the propagation time from spacecraft 2 to 3 and the 
time from 3 to 2: see Fig0] Remark that the maximum 
amplitude of this differential time delay is approximately 

(1/2) 

twice the maximum amplitude of the t contribution, 
because there is a change of sign. Annual variations of 
the frequency shift at order 1/2 are represented on Fig|Sl 
The contribution of order 1 is shown on Fig|H| 



VII. CONCLUSION 

We have given the formulas that allow one to com- 
pute the time delays to be used in a numerical model of 
LISA knowing the Keplerian orbits of the spacecraft. At 
the lowest order, we have the well-known "flexing" of the 
triangle, of (peak to peak) amplitude slightly less than 
48,000 km. At the next order, we find an extra correction 
of amplitude about 960 km, larger than the allowed error 
for phase noise cancellation after TDI. The last computed 



order is negligible (less than 30 m). The interest of the 
derivation is that it includes naturally all relativistic ef- 
fects (for instance the Sagnac effect) introduced up to 
now by "ad hoc" considerations. 

We have also estimated the global frequency shift due 
to the various motions and to propagation in the gravi- 
tational field of the Sun. We have found the rather sur- 
prising result that the General Relativistic contribution 
is much smaller than a priori expectations, only due to 
the LISA configuration. The global frequency shift, re- 
duced to the residual relative motions of the spacecraft, 
is of the order of 8 MHz (peak to peak). 
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FIG. 1: Annual variations of the time delays (zeroth order) along the LISA arms. Solid line: ti2, short dashed line: t23, long 

(0) 

dashed line: tsi. 
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FIG. 4: Annual variations of the Sagnac differential time between spacecraft 2 and spacecraft 3. 



12 



2e-08 




t(s) 



FIG. 5: Annual variations of the contribution of order 1/2 to the frequency shift. 
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